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We study the afterglow phases of a GRB through relativistic magnetohydrodynamic simulations. 
The evolution of a relativistic shell propagating into a homogeneous external medium is followed. 
We focus on the effect of the magnetization of the ejecta on the initial phases of the ejecta-external 
medium interaction. In particular we are studying the condition for the existence of a reverse 
shock into the ejecta, the timescale for the transfer of the energy from the shell to the shocked 
medium and the resulting multiwavelength light curves. To this end, we have developed a novel 
scheme to include non-thermal processeses which is coupled to the relativistic magnetohydrody- 
namic code MRGENESIS in order to compute the non-thermal synchrotron radiation. 
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1. Introduction 



Magnetic fields may play an important role in the relativistic flow of a gamma-ray burst (GRB), 
but the extent to which they are important remains uncertain. Looking at the process which pro- 
duces a relativistic GRB outflow, two alternatives are usually considered. On the one hand, neutrino 
annihilation may be a process which leads to the formation of fireball dominated by thermal energy. 
Here magnetic fields are dynamically unimportant. Alternatively, powerful enough magnetic fields 
can efficiently extract the rotational energy from the central engine and launch a Poynting-flux 
dominated flow. 

In order to achieve relativistic velocities, GRBs have to be launched with high energy-to-mass 
ratio. In the fireball model, most of the energy is assumed to be thermal [17, |^]. This can be be 
result of neutrino-antineutrino annihilation in the polar region of an accreting object [^, |3[ ^, ^. 
The acceleration of the fireball is due to the internal pressure gradient, whereby thermal energy of 
the fireball is converted into kinetic energy of the baryons. At the end of the acceleration phase 
faster parts of the flow may collide with the slower ones producing the internal shocks which power 

|. Processes such as two-stream instability in the shocks 



the GRB prompt emission [31, 11, 27 



might amphfy weak magnetic fields present in the flow []24|]. These fields are expected to account 
for less than 1% of the energy of the flow. After the internal shock phase the flow expands and 
cools as it enters the afterglow phase. In this case very weakly magnetized flow is expected at the 
onset of the afterglow. 

Magnetic fields with appropriate topology can efficiently extract rotational energy from the 
central engine, be it an accretion disk [^, rotating black hole [|^ or a millisecond magnetar [34] 
launching a Poynting-flux dominated flow (PDF). The acceleration of the PDF depends on the 
field geometry and the dissipation processes. Magnetic dissipation can convert Poynting flux into 
kinetic energy [O, 14] and also power the GRB prompt emission [23, 13, 15]. Different studies of 



MHD jet acceleration show that magnetic energy is not completely converted into kinetic energy 
at the end of the acceleration phase. At larger radii we expect magnetic energy of the flow to be 



comparable to the kinetic energy of the baryons [|14j] or even much larger [g3|, It should be 
noted here that the magnetization of the flow might be decreased by the pair-loading caused by the 



vv-annihilation near the central engine [ |22| , 



Fireball and PDF models respectively predict weakly and strongly magnetized flow at the onset 
of the afterglow phase. The initial phases of the interaction of the GRB flow with the (circumburst) 
external medium depend on the strength of the magnetic fields in the flow. A particularly promising 
probe of the magnetization of the GRB flow can come from understanding the early afterglow 



emission [2C, 36, 16]. 



In this paper we outline the status of the ongoing study of the interaction of magnetized ejecta 
with external medium. In Sec. |2| we describe in more detail the current understanding of ejecta- 
medium interaction. Sec. ^ gives an overview of numerical methods and a plan for numerical 
simulations. Summary is given in Sec. H. 
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2. Ejecta-medium interaction 

We consider a homogeneous shell expanding into an external medium of constant density^. 
At large distances from the central engine (typically Rq k, 10^^ — 10^' cm) substantial interaction 
begins, whereby ejecta begins to decelerate due to the accumulation of the external material. We 
assume that at these distances the flow has already been accelerated and coUimated. The internal 
dissipation mechanism, presumably responsible for the prompt emission (e.g. internal shocks, 
magnetic dissipation) is also expected to take place at a shorter distance from the central engine 
with respect to that of the afterglow phases. After the internal dissipation is over, the flow expands 
and cools. Since we are interested in the afterglow phases, the shell is assumed to be cold. We 
denote the shell Lorentz factor by /o ^ 1 and its width by Aq. In a radially expanding outflow 
the magnetic field component perpendicular to the direction of motion drops as r^^ while the 
component in the direction of motion drops as r^^, so that we expect the magnetic field to be 
dominated by the perpendicular component. We define the magnetization parameter as 

where Ep and Ek are Poynting and kinetic energies in the shell, po and Bq its density and magnetic 
field measured in the central engine frame. With this definition a fireball corresponds to ao ^ 1 
while a PDF has ao > 1 . c is the speed of light. The total energy of the shell is 

E = 4nRlAo{YoPoc^ + Bl/4n) = Ek{1 + ao) . (2.2) 



From Eqs. gJ| andg^ we can see that ao parametrizes the fraction of the total energy in the kinetic 
(1/(1 + ao)) and magnetic (ao/(l + ao)) form. 

We first discuss the ejecta-medium interaction for non-magnetized ejecta, and then turn to the 
arbitrarily magnetized case. We also focus on the conditions for the existence of a reverse shock 
into ejecta of arbitrary magnetization. 

2.1 The case of ao< 1 

The evolution of the interface between the cold unmagnetized shell and the external medium 



is well understood. It was studied in detail analytically [32, 26 1 as well as using one-dimensional 



[19] and two-dimensional [ p^ , p^ , 25] numerical simulations. 

At the interface between the shell and the ambient medium two shocks form, the forward 
shock propagating into the external medium, and the reverse shock propagating into the shell. 
Shocked shell and external medium are separated by the contact discontinuity. The forward shock 
is always ultra-relativistic, while the strength of the reverse shock depends on the density contrast 
between the shell and the external medium and the bulk Lorentz factor 70. We distinguish between 



relativistic and Newtonian reverse shocks [|32|]. The critical parameter is 

where / = {3E / AnnempC^y is the Sedov length, ne the external medium number density, and 
rup the proton mass. In the Newtonian case (i§ 3> 1) the shock is non-relativistic in the shell rest 

'similar analysis can be performed for the wind profile where the density of the external medium scales as r^^. 
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frame and does not decelerate the ejecta much, rather the ejecta decelerate once they accumulate 
a mass y^^ times their own mass from the external medium. In the relativistic case (i^ ^ 1) the 
shock crosses the ejecta quickly and slows them down considerably. After the shock crosses the 
ejecta, there is a phase where shocks and rarefaction waves cross the ejecta, passing the energy to 
the forward shock. At later stages the evolution of the ejecta only depends on their total energy and 
the external medium density [0]. 



2.2 The case of arbitrary ao 

The dynamics of magnetized ejecta has not been studied as thoroughly as that of unmagnetized 
ejecta. A qualitative difference to the unmagnetized case is that later phases of the evolution are 
influenced by the internal evolution of the magnetized shell. The initial phase of the evolution has 
recently been studied [ p7[ ] by solving the ideal MHD shock conditions for arbitrarily magnetized 
ejecta with toroidal field. In particular, the dynamics of shock crossing has been studied assuming 
that there is a reverse shock. In that case, the reverse shock crosses the shell faster the higher the 
magnetization is. However, as we have recently shown [[I^, it is not always the case that a reverse 
shock forms. 



2.3 Conditions for the existence of a reverse shock 

Cold, non-magnetized ejecta are always crossed by a reverse shock upon interacting with the 
external medium. This is the case since the sound speed of the ejecta is low and does not allow 
for fast transfer of the information of the interaction with the external medium throughout their 
volume. On the other hand, in a flow that is strongly magnetized and sub-fast magnetosonic (as in 



the Lyutikov & Blandford 2003 model [23]) there is no reverse shock forming. The flow adjusts 



gradually to the changes of the pressure in the contact discontinuity that separates the magnetized 



flow from the shocked external medium. In [16] we generalize to arbitrarily magnetized ejecta and 
derive the condition for the formation of a reverse shock. 

After a detailed treatment of this problem we arrive to the following condition for the formation 



of a reverse shock [16] 



which can be rewritten in terms of shell parameters as^ 

^ A 1/2*3/2 4 T--\I2 ,^ 

ao < 0.6V ^12 72.5^^53 • (2-5) 

Fig. 1 (taken from [|l^]) shows the division of — ao parameter space in two regions, one 
where a reverse shock forms and another where its formation is suppressed. We note that from the 



conditions in Eqs. |2^ and ^ it follows that even for mildly magnetized shells a reverse shock can 
be suppressed. This indicates that the paucity of the observed optical flashes in GRB afterglows 
(associated with the reverse shock emission) may be caused by the suppression of the shock in 
many GRBs. 



^We use the convention that A = A ^ 10' . 
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Figure 1: Existence of a reverse shock in the <^ — (Jq parameter space. The dashed black line delimits 
regions where a reverse shock forms from the region where there is no reverse shock, ignoring the radial 
shell spreading. The solid black line shows the delimitation when the shell spreading is taken into account. 
See [|l6|] for details. 



3. Numerical simulations 

In this section we describe the reasons and motivation for performing numerical simulations 
of shell-ejecta interactions: 



1. Verification of the analytic approach: Results of the analytic work described in Sec. ^, 



especially in Sec. 2.3, need to be verified by means of numerical simulations. We note that 
the line dividing regions of formation and suppression of the reverse shock in Fig. 1, given 



by Eq. 2.4 is approximate, and it is necessary to perform numerical simulations for models 



whose initial parameters lie in the vicinity of the line. 



2. Dynamics of shock propagation: We want to use numerical models to study the influence of 
the magnetization on the propagation of the reverse shock through the shell. On Figs. 2 and 
3 we show snapshots of two simulations, one with the unmagnetized ((Tq = 0), and another 
with the magnetized (ao = 1) shell interacting with the external medium, both shown at the 
same evolutionary time. It can be seen that in the magnetized case (Fig. 3) the reverse shock 
has penetrated the shell deeper than in the unmagnetized case (Fig. 2). Our goal is to perform 
a parametric study where we study the propagation of shocks and rarefactions through the 
shell for different combinations of t, and (Jq. 
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Figure 2: Snapshot from relativistic hydrodynamical simulation of a spherical, non-magnetized (C7 = 0) 
shell that decelerates interacting with external medium of density 10 cm^^. The total energy of the shell is 
E = 10^' erg, its initial width Aq ~ lO'"* cm and bulk Lorentz 70 ~ 50. P* (dashed line), p (solid line) stand 
for the gas pressure and (lab frame) density respectively (in arbitrary units). With increasing radius, one can 
clearly see the reverse shock, contact discontinuity and forward shock located at r ^ 1.483 ■ 10^^ cm. 



3. Long term evolution and energy content: One of the important questions that simulations 
can answer is the timescale of the transfer of energy from the shell to the shocked external 
medium. Long-term numerical calculations are needed to determined the dependence of the 
efficiency of the energy transfer on magnetization. We also want to investigate the long-term 
evolution of the blast wave and determine when the shell profile relaxes to the Blandford- 
McKee solution 

4. Light curves: Finally, we wish to compute synthetic multi-wavelength afterglow light curves. 
Light curves are very sensitive to the magnetic field content of the shell, shock strength and, 
especially, detailed radial profile y{R) of the Lorentz factor as the shell propagates into the 
external medium. They also depend on the distribution of shock accelerated particles. To 
see why y{R) is crucial for the light curve, consider the difference in the arrival time to the 
observer of signals emitted simultaneously from two points with radii R and R + AR, respec- 
tively. It turns out that the difference is At AR y{R)^^ for a relativistic shell. As can be seen, 
a sudden drop in a Lorentz factor in one model as compared to the other will produce features 
which have longer observed duration. The resulting light curves of two models can be very 
different. High-resolution simulations are needed to resolve sufficiently short time intervals 
in order to be able to study the influence of magnetization on y{R) and the subsequent light 
curve. 

We have developed a relativistic magnetohydrodynamic code MRGENESIS [2^ 27, 28] which 
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Figure 3: Same as Fig. 2, but for the magnetized (do ~ 1) shell. Note that the reverse shock crosses the 
ejecta faster with respect to the non-magnetized case in agreement with analytical expectations. 



consists of a finite- volume, high-resolution, shock-capturing scheme GENESIS @, 21 1 which 
solves for the conservation laws of relativistic magnetohydrodynamics, and a module which fol- 
lows the transport, evolution and radiation from non-thermal particles. For the purpose of analyzing 
dynamics and emission from afterglow shells, a very high resolution is needed. Numerical conver- 
gence tests have shown that, in order to sufficiently resolve the shell, we need to resolve the scales 
of the order of Ao7o"^. ^^^^ means that we need to use at least zones within the shell. We use a 
grid re-mapping procedure described in |29], which can be thought of as a guided mesh refinement 
centered on the shell. Even with this procedure the actual number of zones for a typical model is 
expected to be several millions. We are performing simulations on supercomputers of the Spanish 
Supercomputing Network. 



4. Conclusions 

We are performing high-resolution numerical studies of the transition from prompt to the early 
afterglow phase of gamma-ray bursts. The afterglow is being modeled as the radiation from the rel- 
ativistic shell expanding into the homogeneous external medium. We study the difference between 
the fireball (unmagnetized shell) and the Poynting-flux dominated (magnetized shell) models. 

In the context of the early optical afterglow we have shown analytically that even a moder- 
ate magnetization of the flow can suppress the existence of a reverse shock, and thus explain the 
apparent paucity of the optical flashes for a large number of early afterglows. We are currently 
performing simulations to study the formation and suppression of relativistic shocks in detail. To 
study the later phases of the afterglow we aim to determine the influence of the initial shell mag- 
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netization on the energy content, transfer of energy from the shell to the forward shock, and the 
long-term flow structure. 

We have developed a novel scheme for treating non-thermal processes in relativistic magne- 
tohydrodynamic simulations. This scheme is used to compute multi-wavelength light curves from 
numerical simulations. The aim is to study the influence of the initial magnetization on the short- 
and long-term Ught curves. 
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